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^ , , We study the exclusive decays of B meson into pseudoscalar charmonium states 

O ' r]c and r]'^ within the QCD factorization approach and find that the nonfactorizable 

■ corrections to naive factorization are infrared safe at leading-twist order. The spectator 

\ interactions arising from the kaon twist-3 effects are formally power-suppressed but 

I chirally and logarithmically enhanced. The theoretical decay rates are too small to 

^ ■ accommodate the experimental data. On the other hand, we compare the theoretical 

. calculations for J/'0,V'', and r]c,r]'^, and find that the predicted relative decay rates of 

I these four states are approximately compatible with experimental data. 

o\ : 

O . PACS numbers: 13.25.Hw; 12.38.Bx; 14.40. Gx 

(N ■ 
O ; 

1 Introduction 

Exclusive decays of B meson to charmonium are important since those decays e.g. B — >■ 
^ ■ J/ ipK are regarded as the golden channels for the study of CP violation in B decays. How- 
>• . ever, quantitative understanding of these decays is difficult due to the strong-interaction 
^ ! effects. It is conjectured physically that because the size of the charmonium is small(~ 
^ I l/agm^) and its overlap with the {B,K) system is negligible the same QCD-improved 
factorization method as for i? ^ 7r7r[0, ^ can be used for B —>■ J/ipK decay. Indeed, for 
this channel the explicit calculations 0, Q show that the nonfactorizable vertex contribution 
is infrared safe and the spectator contribution is perturbatively calculable at twist-2 order. 
This small size argument for the applicability of QCD factorization for the charmonia is 
intuitive, but it needs verifying for charmonium states other than the J/ijj and ip'{ip{2S)). 

In our previous paper 0, we studied the B — > XcjK{J = 0, 1) decays within the QCD fac- 
torization approach, and found that for B — > XaK decay, the factorization breaks down due 
to logarithmic divergences arising from nonfactorizable spectator interactions even at twist-2 
order, and that for B XcqK decay, there are infrared divergences arising from nonfactoriz- 
able vertex corrections as well as logarithmic divergences due to spectator interactions even 
at leading- twist order. 

Experimentally, for the pseudoscalar charmonium state r^c? the B rj^K decay has 
been observed by CLEOjj^, BaBar|p, and Belle ||^ with relatively large branching fractions. 
Moreover, very recently the ri'Jj](.{2S)) meson has also been observed in the B —>■ rj'^K decay 
by Belle|TDf|. So, it is interesting to compare the predictions of these decay modes into 
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pseudoscalar charmonium based on the QCD factorization approach with the experimental 
data to further test the apphcabihty of QCD factorization to B meson exclusive decays to 
charmonium states. 



2 B ^ 7]cK decay within QCD factorization 

We now consider B rjcK decay. The effective Hamiltonian for this decay mode is written 



as 
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V / i=3 

where Ci are the Wilson coefficients and the relevant operators Oi in Hcs are given by 

Oi = {Sabp)v-A ■ {Ci3Ca)v-A, O2 = (s„6„)v-A " {Cf3Cp)v-A, 

^3(5) = {Saba)v-A ' Y{Qi3(1i3)v^A{V+A), (^4(6) = {Sabi3)v-A " Y{Qi3(la)v^A{V+A), (2) 

3 3 

^7(9) = 7^i'^aba)v-A " 'Yeq{Qi3Ql3)v+A{V-A), C'sCiO) = 7^ i'^ ab /s) V - A " eq{Qi3la)v+AiV-A)- 
^ q ^ q 

To calculate the decay amplitude, we introduce the rjc decay constant as 

(^c(p)|c(0)7;.75c(0)|0) = -tUp^., (3) 

where /^^ is the rjc decay constant which can be estimated from the QCD sum rules or 
potential models. The leading-twist light-cone distribution amplitude of rjc is then expressed 
compactly as 

{Vcip)\caiz2)cpizm = ^j^'c^«-e^("^-^+(i-")^-^^)[(|^ + m,j75],>,.H, (4) 

where u and 1 — u are respectively the momentum fractions of the c and c quarks inside the 
rjc meson, and the wave function for 7]c meson is symmetric under u ^ 1 — u. 

As for the kaon light-cone distribution amplitudes, we will follow Ref. Qto choose 

{K{p)\s^{z2)d^iz,)\0) 
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where x and 1 — x are respectively the momentum fractions of the s and d quarks inside 
the K meson. The asymptotic limit of the leading-twist distribution amplitude is (pKix) = 
6x(l — x). We also use the asymptotic forms = 1 and = 6x(l — x) for the 

kaon twist-3 two-particle distribution amplitudes. The chirally-enhanced factor is written 
as = 2fiK/mh = 2mj^/mb{ms + md) which is formally of order Aqco/^^fe but numerically 
close to unity. 

In the naive factorization, we neglect the strong interaction corrections and the power 
corrections in Aqcd/'^6- Then the decay amplitude is written as 

^Mo = ^UmlFo{ml)^[v^,V:,{C, + ^) - V,,V,:{C, + ^^-C,-^)], (6) 
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Figure 1: Feynman diagrams for nonfactorizable corrections io B rjcK decay. 



where Nc is the number of colors. We do not include the effects of the electroweak penguin 
operators since they are numerically small. The form factors for B K are given as 



{K{pK)\sl^J:)\B{pB)) = {Vb+Pk, 
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where p = Pb ~ Pk is the momentum of r]c with = m^^ and we will neglect the kaon mass 
for simplicity. 

As we can see easily in Eq. this amplitude is unphysical because the Wilson coeffi- 
cients depend on the renormalization scale /i while the decay constant and the form factors 
are independent of /i. This is the well known problem with the naive factorization. How- 
ever, if we include the order corrections, it turns out that the /i dependence of the Wilson 
coefficients is cancelled and the overall amplitude is insensitive to the renormalization scale. 
Taking the nonfactorizable order strong- interaction corrections in Fig. |I| into account, the 
full decay amplitude for B rj^K within the QCD factorization approach is written as 

G p 

where the coefficients a, {i = 2, 3, 5) in the naive dimension regularization(NDR) scheme are 
given by 
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Ci 


C2 


C3 


C4 


C5 


Ce 


LO 
NDR 


1.144 
1.082 


-0.308 
-0.185 


0.014 
0.014 


-0.030 
-0.035 


0.009 
0.009 


-0.038 
-0.041 



Table 1: Leading-order(LO) and Next-to-leading-order(NLO) Wilson coefficients in NDR 
scheme(See Ref.Q) with = 4.4 GeV and Ag^ = 225 MeV. 
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The function fj in Eq.(|^) is calculated from the four vertex corrections (a,b,c,d) in Fig. |T| 
and it reads 
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where z = m,^^/m^, and we have already symmetrized the result with respect to 

The function /// in Eq.(P) is calculated from the two spectator interaction diagrams (e,f) 
in Fig. and it is given by 



B 



Nr m%Fn(ml 







du 



(prjciu) dx 



U 



— [(PK{X} 
X 



2nK<pKix) , 

mb(l - z) 



(11) 



where (ps is the light-cone wave functions for the B meson. The spectator contribution 
depends on the wave function (ps through the integral 
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Since (psiO is appreciable only for ^ of order Aqcd/'^b, is of order Aqcd- We will follow 
Ref. [P] to choose ~ 300 MeV in the numerical calculation. 

There is an integral in Eq. (^) arising from kaon twist-3 effects, which will give logarith- 
mic divergence. Following Ref. 0, we treat the divergent integral as an unknown parameter 
and write 

"1 dx 
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where 0^(a;) = 1 is used for the kaon twist-3 light-cone distribution amplitude. We will 
choose Xh = iTL^niB/ Aqcd) ~ 2.4 as a rough estimate in our calculation. 



For numerical analysis, we choose Fo(m^J 



0.41 [|T2] and use the following input param- 



eters: 



TUb = 4.8 GeV, niB = 5.28 GeV, m„^ = 3.0 GeV 



fr,, = 350 MeVim, fs = 180 MeV, fx = 160 MeV. 



(14) 
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[u) 


Ci2 






as 


Qu{l 


-u) 


0.1043-0.06^ 


Mi 


0.0045+0.0022i 


-0.0035-0.0026i 


6{u~ 


1/2) 


0.0792-0.06^ 


I2i 


0.0055+0.0022i 


-0.0048-0.0026i 



Table 2: The coefficients a^ at /x = 4.4 GeV with different choices of 



The asymptotic form of the distribution amphtude is given as = 6-u(l — -u). 

In the numerical analysis, we also consider the form (prj^iu) = 5{u — 1/2), which comes 
from the naive expectation of the distribution amplitude. Although there are uncertainties 
associated with the form of the wave function, we will see shortly that the calculated decay 
rates are not very sensitive to the choice of the distribution amplitude. The results of 
coefficients a^ are listed in Table. ^ 

With the help of these coefficients a^, we calculated the decay branching ratios. For 
<pr,Xu) = Qu{l - u), Bi(B r]^K) = 1.9 x 10^1 And for (f)^^^{u) = 6{u - 1/2), Bi(B 
r]cK) = 1.4 X 10-^ 

The measured branching ratios are 

CLEO Collaboration 0: Br(fi° ^ r]^K^) = (1.091°;^^) x 10-^ 
BaBar Collaboration §: Br(5° r]^K^) = (1.06 ± 0.28) x 10-^ (15) 
Belle Collaboration 0: Br(5° r]cK^) = (1.23 ± 0.23) x 10~^ 

which are about seven times larger than our theoretical results. 



3 B ^ q'^K decay 



The calculation of the branching ratio for B rj'^K decay is similar to that for the rjc given 
above. And we can also get a rough estimate for the decay rates ratio of rj'^ to rjc in the 
leading order: 

Br(B"^r?c^") /^^ Fxim^^J m% - m^^^ /^^ 

where we have used fn'Jfri^ ~ fip'/ fj/-^ with /j/^, = 400 MeV, = 280 MeV, which are 
determeined from the observed leptonic decay widths and -Po(p^) = (1 — / fTi%)Fi{p'^) 



with Fi{m\) = 0.81,Fi(m2j = 0.58 [|, The ratio in Eq.(|l|) will roughly hold even 
when we include the 0{as) corrections, because the 0{as) corrections are small and the 
mass difference as well as the wave function difference between r]c and rj'^ will not change the 
values of a^ in Eq.(|^) greatly. 

The Belle Collaboration has reported the observation of the r]'^ in exclusive B —>■ KKsK^tt^ 
decays ||T0[]: 

Br(5° ^ 7]'^K^)Bi{7]'^ KsK-TT+ 



Br(50 ^ r/cisTO) Br (r/c KsR-n^ 



0.38 ±0.12 ±0.05. (17) 
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As was noted in Ref.||15||, the hadronic decay branching fractions for t]c and 77^ are expected 



to be roughly equal for the helicity non-suppressed decay channels^. So we have Br(?7^ 
KsK~7i~^) ^ Br(?7c —>■ KsK~ 7r~^), and then from Eq.([T7|) we get 



^ lkLLi~04 CIS) 

which is consistent with the ratio in Eq.(|IB|). However, as was mentioned above, because 
the theoretical decay rate is about seven times smaller than the the experimental data for 
Br(i?° — > rjcK^), the theoretical branching fraction will also be about seven times smaller 
than the the experimental data for Br (5° —>■ rj'^K^). 



4 Discussion 

We have shown that for B decays to rjc and 1]'^ the theoretical branching fractions are all about 
seven times smaller than the experimental data. However, from Eq.(|T^) and Eq.(|TBp, we see 
that the theoretical ratio of the decay rates of the two states is consistent with experimental 
data: 

Br(i?° ^ i^K^) ^ Br(i?° ^ 

Br(50^r/,iro)^^/^ Br(50 ^ r/,iro) ' ^ ' 

It is also interesting to find that although the theoretical branching fractions of B meson 
exclusive decays to J/ and tp' are both much smaller than the experimental data, the theo- 
retical ratio of the decay rates of these two states is also roughly consistent with experimental 
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V I ^ ( jt^ ^2 ^ i\ I ^2 ^ B i, i3 _ Q 9 X (_^)2 ~ 0.45, (20) 



Br(5° ^ i)'K^) 

MB^^J/i^K^)^^ ^ ^^^^ 

where Fiim^,) = 0.83 and Fi(mj^^) = 0.61 are used. 

Another interesting observation is that the theoretical ratio of the branching fractions 
of B meson exclusive decays to rjc and J/ip is also roughly consistent with experimental 
data0, 1, 1, IH: 



V /c J f J 'i'^' _ I ^ Vc' ]^ I o ^'^i~10x i ) ~ 75 (22) 



Br(50 ^ r7cK°) , , 

Br(50^ J/V^i^o)^^ ■ ■ ^ ^ 



^It will also be interesting to detect the helicity suppressed decay channels of rjc and ry^ in B decays, 
and to see the differences between the helicity suppressed (e.g. pp, K* K* , cfxp pp) and non-suppressed (e.g. 
KKTT,ri7nr,ri'7nT) decays of rjc and 77^. This will be useful to clarify the helicity suppression mechanism for 
the charmonium hadronic decays and the so-called pn puzzle in J/i p and ip' decays observed by BES and 
annihilation experiments. For details, see Refs.||15|, |l6|. 



6 



Eq.(p2D also approximately holds when 0{as) corrections are included. 

So, the predicted relative rates of all S-wave charmonium states J/ip, r]c, Vc QCD 
factorization approach are roughly compatible with data. This has been shown explicitly 
above in the leading order approximation, and even holds when including 0{as) corrections 
with which the calculated decay rates for these four charmonium states are almost equally 
smaller than data by a factor of 7-10 though there are some theoretical uncertainties associ- 
ated with form factor^, decay constants, as well as the light-cone wave functions of mesons 
involved. This result is rather puzzling, and it might imply that the naive factorization for 
B decays to the 5'-wave charmonia may still make sense but the overall normalization for 
the decay rates are questionable. 

In summary, we have studied the exclusive decays of B meson into pseudoscalar charmo- 
nium states 7]c and r]'^ within the QCD factorization approach and find that the nonfactoriz- 
able corrections to naive factorization are infrared safe at leading-twist order. The spectator 
interactions arising from the kaon twist-3 effects are formally power-suppressed but chirally 
and logarithmically enhanced. The theoretical decay rates are too small to accommodate the 
experimental data. We already knew that for B J/ipK decay, there are also logarithmic 
divergences arising from spectator interactions due to kaon twist-3 effects and the calculated 
rates are also smaller than data by a factor of 8-100, ^. Moreover, in our previous paper 
we found that for B XciK decay, the factorization breaks down due to logarithmic di- 
vergences arising from nonfactorizable spectator interactions even at twist-2 order, and the 
decay rates are also too small to accommodate the data, and that for B XcoK decay, there 
are infrared divergences arising from nonfactorizable vertex corrections as well as logarithmic 
divergences due to spectator interactions even at leading-twist order. 

Considering the above problems encountered in describing B —>■ rjcKirj'^K) as well as -B — >■ 
J /ipK{ip' K) and B XcjK{J = 0, 1) decays, we would like to restate our conclusion that 
in general the QCD factorization method with its present version can not be safely applied 
to exclusive decays of B meson into charmonia, and that new ingredients or mechanisms 
should be introduced to describe exclusive decays of B meson to charmonium states. 



Acknowledgements 

We thank G. Bodwin, E. Braaten, X. Ji and C.S. Lam for useful discussions and comments. 
This work was supported in part by the National Natural Science Foundation of China, and 
the Education Ministry of China. 

^In our calculation, we have used the relation -Fo(p^) — (1 ~ p'^ /'m%)Fi{p'^) derived in Ref.[^, which is 
consistent with the form factors obtained in Ref . [p^ . This will reduce the effects of uncertainties arising 
from form factors on the decay rate ratios. For example, in Eq.(^2|) we have used 

(1 - ml/minp^r « (1 - ml/mir = 0.46. 

This value is close to that given in Ref. , which is the modified version of Ref. ||l^ . In Ref. , the authors 
discussed the B decay rate ratio of r]c to J/ip at the leading order and assumed that -Fo(p^) is a constant 
and Fi{p'^) has a monopole dependence with specific pole masses. 
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